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Abstract 

The evolution of neutron stars in close binary systems with a low-mass companion is 
considered assuming the magnetic field to be confined within the solid crust. We adopt the 
standard scenario of the evolution in a close binary system in accordance with which the 
neutron star passes throughout four evolutionary phases ("isolated pulsar" - "propeller" - 
accretion from the wind of a companion - accretion due to Roche-lobe overflow) . Calculations 
have been performed for a great variety of parameters characterizing the properties both of 
the neutron star and low-mass companion. We find that neutron stars with more or less 
standard magnetic field and spin period being processed in low-mass binaries can evolve to 
low-field rapidly rotating pulsars. Even if the main-sequence life of a companion is as long 
as 10 10 yr, the neutron star can maintain a relatively strong magnetic field to the end of the 
accretion phase. The considered model can well account for the origin of millisecond pulsars. 

Key words: pulsars: recycled - neutron stars: structure - neutron stars: magnetic fields 
- X-ray binaries 

1 Introduction 



In recent years evidences have been obtained for a relatively fast magnetic field decay in neu- 
tron stars undergoing accretion in binary systems (see, e.g., Taam & van den Heuvel 1986, 
Bhattacharya & van den Heuvel 1991). The weak surface magnetic fields of many binary ra- 
diopulsars (B < 10 9 G), for which an efficient accretion has taken place, strongly suggest the idea 
of the accretion induced magnetic field decay in these neutron stars. Most of the low- magnetic- 
field pulsars are believed to be old neutron stars experienced a mass transfer in low-mass binaries 
where accretion ensues due to Roche-lobe overflow and may last relatively long. 

Recently Geppert & Urpin (1994) and Urpin & Geppert (1995) have suggested a simple 
and efficient mechanism of a rapid field decay in neutron stars undergoing accretion. Accretion 
changes drastically the thermal evolution of the neutron star if the accretion rate is sufficiently 
high and a mass transfer lasts a sufficiently long time (see, e.g., Fujimoto et al. 1984). The 
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energy release due to accretion heats the neutron star and reduces the crustal conductivity since 
the latter depends on the temperature. A decrease of the conductivity accelerates the ohmic 
decay of the magnetic field if this field is maintained by currents in the crust. If the mass 
transfer phase lasts as long as 10 6 — 10 7 yr the surface field strength can be reduced by a factor 
~ 10 2 — 10 4 depending on the magnetic configuration at the beginning of the accretion phase. 

Accretion influences not only the field decay but the spin evolution of pulsars as well, and the 
strength of the magnetic field is one of the dominating factors that determines this evolution. 
If the total amount of the accreted matter is large enough (> 0.01 — 0.1M o ), the neutron star 
in a low- mass binary can be "recycled" to millisecond periods (Alpar et al. 1982). 

The accretion phase, however, does not exhaust all possible stages in the evolution of neutron 
stars entering low-mass binary systems. The evolution of such systems is extremely long since 
even the main-sequence lifetime of a low-mass companion exceeds 10 9 yr. In accordance with the 
standard scenario, the neutron star in a close binary can pass throughout several evolutionary 
phases (see, e.g., Pringle &: Rees 1972, Illarionov & Sunyaev 1975): 

1) the initial obscured radio pulsar phase in which the pressure of the pulsar radiation is 
sufficient to keep plasma of the companion stellar wind away from the neutron star magneto- 
sphere; the radio emission of the pulsar is partly (or completely) absorbed by the plasma cloud 
surrounding the binary system thus the pulsar is practically unobservable; the pulsars spins 
down due to the magnetodipole radiation; 

2) the propeller phase in which the radiation pressure reduced due to the field decay and 
spin down cannot prevent the infalling matter from an interaction with the magnetosphere but 
the spin rotation is sufficiently rapid yet to eject the wind matter due to the centrifugal force 
and to transfer the angular momentum from the neutron star to the wind matter; 

3) the wind accretion phase in which the matter from the wind falls on to the surface of the 
neutron star; nuclear burning of the accreted material heats the neutron star and causes the 
accretion-driven field decay; in the course of this phase the neutron star can slightly spin up 
since the accreted matter carries a certain amount of the angular momentum; 

4) the enhanced accretion phase which starts when the secondary star leaves the main- 
sequence and fills its Roche-lobe; a heavy mass transfer on to the neutron star heats the interior 
to a very high temperature ~ 10 8 — 10 8 5 K (Fujimoto et al. 1984, Miralda-Escude et al. 1990) and 
accelerates the accretion-driven field decay; a steady Keplerian disc is formed by the accretion 
flow outside the neutron star magnetosphere and the accretion torque spins up the neutron star 
to a very short period. 

The neutron star processed in the above evolutionary picture can probably prolong its life 
as a radiopulsar with relatively low magnetic field and short spin period after accretion due to 
Roche-lobe overflow is exhausted. The properties of pulsars produced due to evolution in low- 
mass binaries have not been analysed in detail yet and we address this problem in the present 
paper. Note that a similar problem has been considered by Jahan Miri & Bhattacharia (1994) 
who assumed that the magnetic evolution is determined by the expulsion of the field lines from 
the neutron star core due to interaction between the fluxoids in the proton superconductor and 
the vortices in the neutron superfluid. The authors, however, restricted themselves to the case 
of wide low-mass binaries and considered the evolution only during the main-sequence life time. 

The paper is organized as follows. We describe in detail the adopted model in Section 2. 
The results of calculations are presented in Section 3. In Section 4, we discuss our results and 
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compare them with the available observational data on the evolution of neutron stars in low-mass 
binaries. 



2 Description of the model 

We consider the evolution of the neutron star in a binary system with the secondary being a 
low-mass star. The binary is assumed to be relatively close thus the companion can fill the Roche 
lobe after the end of its main-sequence life. Both the magnetic and spin evolution of the neutron 
star may be strongly influenced by accretion from the companion. During its main-sequence 
evolution (which lasts as long as (1 — 5) x 10 9 yr), the secondary loses mass due to the stellar 
wind and a certain fraction of the wind plasma can generally be captured and accreted by the 
neutron star. After the secondary leaves the main-sequence and fills the Roche-lobe, the rate of 
mass loss and, hence, the accretion rate can be drastically enhanced. 

We assume that the magnetic field of the newly born neutron star is created by an unspecified 
mechanism either in the course of collapse or soon after the neutron star is born. We also assume 
that the field is maintained by currents in the crust and occupies originally only a fraction of the 
crust volume. The evolution of the crustal field is determined by the conductive properties of 
the crust and material motion throughout it (see for a detail Sang & Chanmugam 1987, Urpin 
&; Geppert 1995). Neglecting the anisotropy of conductivity, the induction equation governing 
the evolution of the magnetic field reads 



where a is the conductivity and v is the velocity of material movement. The velocity v is caused 
by the flux of the accreted matter throughout the crust and may be non-zero only during those 
evolutionary phases when accretion on to the neutron star surface is allowed. The total amount 
of the accreted mass in most cases does not exceed 0.1-0. 2M (see Bhattacharya & van den 
Heuvel 1991), even for neutron stars undergoing a very extended accretion phase. Therefore, 
the mass and radius of the neutron star do not change significantly and, hence, the corresponding 
changes in the crustal thickness, density profile etc. are also negligible. The total mass flux is the 
same throughout the crust and should be taken equal to the accretion rate. Assuming spherical 
symmetry of the mass flow in deep crustal layers, one obtains the radial velocity of material 
movement 



where M is the accretion rate and p is the density. This inward directed flow tries to push the 
magnetic field into the deep crustal layers and, generally, can influence the evolution of the field 
at a high accretion rate. 

Introducing the vector potential for a dipole magnetic field A = (0,0, A v ), where A v = 
s(r,t) sin 9/r and r,9,ip are the spherical coordinates, the induction equation (1) can be trans- 
formed to 




(1) 




(3) 
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At the surface r = R, the standard boundary condition for a dipole field should be fulfilled, 
i.e. Rds/dr + s = 0. We assume the core of the neutron star to be superconductive, thus the 
magnetic field cannot penetrate into the core. Therefore, the second boundary condition is s = 
at the crust-core boundary which lies at the density 2 x 10 14 g/cm 3 . 

The conductive properties of the solid crust are determined by scattering of electrons on 
phonons and impurities (see, e.g., Yakovlev & Urpin 1980). Scattering on phonons dominates 
the transport processes at a relatively high temperature and low densities, whereas scattering 
on impurities gives the main contribution to the conductivity at a low temperature and large 
densities. The total conductivity of the crystallized region is given by 



where a p h and <r imp are the conductivities due to electron-phonon scattering and electron- 
impurity scattering mechanisms, respectively. The phonon conductivity depends on the tem- 
perature, decreasing when T increases, and hence heating caused by accretion accelerates the 
field decay. The impurity conductivity is practically independent of the temperature but its 
magnitude is determined by the so called impurity parameter, 



where the dominant background ion species has the number density rii and charge Z, and n' 
is the number density of an impurity species of charge Z'\ the summation is over all species of 
impurities. In our calculations, we use the numerical data for the phonon condutivity obtained 
by Itoh et al. (1993) and a simple analytical expression for the impurity conductivity derived 
by Yakovlev & Urpin (1980). Since the conductivity depends on the temperature, the magnetic 
evolution of the neutron star is strongly sensitive to its thermal history. According to the sce- 
nario suggested in Section I, we divide the evolution into four essentially different phases: 
the "isolated pulsar" phase (phase I) in which the neutron star practically does not feel the 
influence of its companion, 

the "propeller" phase (phase II) in which, due to a fast spin rotation, the neutron star ejects 
the wind matter penetrating into a magnetosphere, 

the wind accretion phase (phase III) in which accretion of the wind plasma is allowed, and 
the accretion phase (phase IV) in which the secondary fills its Roche-lobe and the neutron star 
experiences a heavy mass transfer. 

Evidently, that during the phases I and II the thermal evolution of the neutron star in a 
binary system exactly follows that of an isolated neutron star. For the time dependence of the 
temperature during these two phases we use the standard cooling model as given by Van Riper 
(1991) for 1AM & neutron star with the Pandharipande-Smith equation of state. During the 
phases III and IV, however, the neutron star can be substantially heated by accretion. The 
thermal evolution of accreting neutron star has been a subject of study of several papers (see, 
e.g., Fujimoto et al. 1984, Miralda-Escude et al. 1991, Zdunik et al. 1992). It was argued 
that, after a relatively short initial stage, the temperature distribution within the neutron star 
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reaches a steady state. Then, the temperature is practically uniform in deep layers of the crust 
of neutron star with standard neutrino emissivities. The temperature that the crust will attain 
in the steady state depends on the accretion rate, M. For the phase III with a low accretion 
rate, we use the crustal temperature computed by Zdunik et al. (1992) for the accretion rate 
within the range 1CP 11 > M > lO _16 M0/yr. For the phase IV when accretion is caused by 
Roche- lobe overflow, the internal temperature has been calculated by Fujimoto et al. (1984) 
for the accretion rate M = 3 x 10~ 10 M Q /yr and by Miralda-Escude et al. (1991) for M within 
the range 1CT 10 > M > lO~ 11 M /yr. We matched their results and slightly extrapolate the 
obtained dependence in order to estimate T for M ~ lO -9 M0/yr, the maximal accretion rate 
we use in calculations. 

Note that the accretion rate during the phases III and IV is sensitive to the separation 
between the stars and their masses. In the present study we will not consider the dependence of 
the evolution on the orbital period of a binary despite of the importance of this dependence for 
interpretation of observational data (see van den Heuvel & Bitzaraki 1995). We are planning to 
address this problem in a forthcoming paper. 

The spin evolution of the neutron star being strongly coupled with its magnetic evolution 
differs essentially for different evolutionary phases. As it was mentioned, the neutron star is 
practically not influenced by a companion during the phase I since the wind plasma of the 
secondary is stopped by the pressure of magnetodipole radiation behind the light cylinder. The 
stopping radius, R s , is determined by the balance between the dynamical pressure of the wind 
and the radiative pressure of magnetodipole waves. The dynamical pressure of the wind is of the 
order of p w Vw where p w and V w are the density and velocity of the wind plasma, respectively, 
near the stopping radius. The density, p w , at the distance s from the secondary may be estimated 
as p w f» Mq/4,tts 2 V w with Mq being the rate of mass loss of the secondary. The pressure of the 
magnetodipole radiation at the distance R s from the neutron star is ~ (l/47rR 2 )(B 2 R G ft 4: /6c 4: ) 
where B is the surface field strength at the magnetic pole and £1 is the spin angular velocity of 
the neutron star. Equating this pressure to the dynamical pressure of the wind at the nearest to 
the neutron star point and assuming that the stopping radius is smaller than separation between 
stars, a, i.e. s ~ a, one obtaines the expression for R s , 

„ (b'r'W o 2 \ 1/2 

In our calculations, it is more convenient to use the accretion rate, M, instead of the rate of 
mass loss Mq. At a > Rg where Rg ~ 2GM/V 2 is the radius of gravitational capture (Bondi 
1952) and M is the neutron star mass, these quantities are related by 4M/R G « Mo/a 2 ; we 
assume that the wind velocity is larger than both the sound velocity of the wind plasma and the 
velocity of orbital motion. The phase I lasts while the stopping radius is larger than the radius 
of gravitational capture, R s > Rg- 

With the assumption that the spin down torque on the neutron star is caused by radiation 
of the magnetodipole waves (Ostriker & Gunn 1969), the evolution of the spin period, P, is 
governed by 

• 8.r 2 B 2 iJ 6 
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where / is the moment of inertia. Note that approximately the same braking law applies if the 
torque is determined by currents in the model suggested by Goldreich and Julian (1969). At 
R s < Rg, the pressure of the magnetodipole radiation cannot prevent the wind plasma from 
a penetration behind the light cylinder. The influence of the wind plasma on the spin-down 
rate becomes appreciable and the further spin evolution of the neutron star departs from that 
of an isolated pulsar. In our simplified model, the condition R s = Rq determines the transition 
between the phases I and II. 

During the phase II, the wind matter can directly interact with the neutron star magne- 
tosphere. However, rotation is rather fast yet and the magnetosphere acts as a "propeller", 
ejecting the infalling wind matter. During this phase, a certain fraction of the spin momentum 
of the neutron star is transfered to the matter being ejected, and the spin period may increase a 
few orders of magnitude. In calculations presented here, we adopt a very simplified model of the 
interaction between the accretion flow and magnetosphere which is, nevertheless, qualitatively 
close to that commonly encountered in the literature (see, e.g., Pringle & Rees 1972, Illarionov 
& Sunyaev 1975). We assume that the accreted matter interacts with the magnetosphere at the 
so called Alfven radius, Ra, 

( r6 R 2 \ 2 / 7 / R 2 \ 2 / 7 

' ' ■ < « L2x i 9 cm, (8) 



A \4MVGm) 



where B\2 = -B/10 12 G and M_io = M/lO~ lo M0yr _1 ; the numbers are given for the neutron 
star with M = 1AM & and R = 1.64 x 10 6 cm. If the neutron star rotates rapidly and its angular 
velocity is larger than the Keplerian angular velocity at the Alfven radius, Qk(Ra), the wind 
matter has to be expelled by the magnetosphere carrying out the spin angular momentum. In 
order to be expelled, the wind matter should get the angular momentum a little bit larger than 
it corresponds to the Keplerian rotation. Therefore, the amount of the spin momentum lost by 
the neutron star per unit time, J p , may be estimated as 

j p ~ 4-irR\p w n K v r ~ M^GMR A , (9) 

where v r is the radial velocity of the accreting flow near the magnetospheric boundary. The 
corresponding spin-down rate of the neutron star is 

p~El±~ f3P 2 B 2 l 7 M & l 7 , (10) 
2irl 

where (5 = (GMR 2 /8) 3 ^ 7 /irl. The neutron star rotation slows down until the angular velocity 
becomes comparable to Uk(Ra)- We determine the critical period, P eq , at which the transition 
between the "propeller" and accretion phases takes place from the condition Q = Qk{Ra), 

ivi_ w 

This period can be rather large if the neutron star is accompanied by a secondary star with a 
low rate of mass loss. Equation (11) determines the so called spin up line in the B - P plane. 
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If the spin rotation slows down to the value given by equation (11), it is generally accepted 
that the accreting matter can flow throughout the magnetosphere and accrete on to the surface, 
thus the neutron star enters the phase III. During this phase accretion from the stellar wind 
is allowed and, as it was mentioned, nuclear burning of the accreted material provides some 
heating of the neutron star interior and induces the accretion-driven field decay. In its turn, the 
decay of the magnetic field decreases the critical period, P eq . Since accretion from the wind in 
a binary cannot be exactly spherical, the accreting matter carries some amount of the angular 
momentum to the neutron star (see Geppert et al. 1995) and, due to this, the neutron star can 
experience spin up to a shorter period. Probably, a balance may be reached in spin up and the 
rate of the field decay, thus the neutron star slides down the corresponding spin-up line (see, e.g., 
Bhattacharya & Srinivasan 1991). This evolution lasts either until the end of the main-sequence 
lifetime of the companion (when it fills the Roche-lobe and the neutron star enters the phase IV) 
or until the magnetic field becomes weak to such extent that the accreting angular momentum 
is insufficient to maintain a balance in spin up and the rate of the field decay. 

The amount of the angular momentum carried by the accreted matter can be estimated from 
the condition at the magnetospheric boundary assuming that the angular momentum there is 
characterized by its Keplerian value, pVLkR\, multiplied by a some "efficiency" factor £. This 
factor is ~ 1 if the accreted wind matter forms a Keplerian disc outside the magnetosphere. 
If the disc is not formed, then the efficiency of spin up is much lower and the parameter £ is 
probably of the order of 0.1 - 0.01 (see Jahan Miri & Bhattacharya 1994). The amount of the 
spin angular momentum added to the neutron star per unit time, J w , can be estimated as 

j w ~ £ • 4irR 4 AP n K v r ~ £M^GMR A . (12) 

The corresponding spin-up rate is 

P ~ -£PP 2 B 2/7 M 6/7 . (13) 

If the magnetic field becomes very weak, P may be small and accretion cannot maintain a 
balance in the spin up and field decay rates, thus the neutron star leaves the spin-up line. 

After the companion fills its Roche-lobe and the phase IV starts, accretion onto the neutron 
star is strongly enhanced and heating of the crust leads to a rapid field decay. In low-mass 
binaries, accretion due to Roche- lobe overflow can probably last as long as (1 — 5) x 10 7 yr. 
During this phase, the accreted matter likely forms the Keplerian disc outside the neutron star 
magnetosphere. Therefore, the rate of the angular momentum transfer and the spin-up rate 
are given by equations (12) and (13), respectively, with the "efficiency" parameter £ ~ 1. The 
neutron star spins up in accordance with these equations until it approaches the spin-up line 
corresponding the enhanced accretion rate. Probably, the accretion rate has to reduce when the 
neutron star reaches the spin-up line (Bhattacharya & Srinivasan 1991) but, nevertheless, the 
field continues to decay due to the accretion-driven mechanism. During the further evolution, 
a balance should be reached in spin up and the rate of the field decay like the case of accretion 
from the wind. Due to this balance, the neutron star slides down the corresponding spin-up line 
with the rate which is determined by the field decay. In reality, the evolution on the spin-up 
line may be very complicated but, in our calculations, we will adopt a maximally simplified 
model. Namely, we will neglect a possible difference in the accretion rate before and after the 
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neutron star reaches the spin-up line. As mentioned above, the amount of matter accreted onto 
the surface may be reduced when the star evolves on the spin-up line. However, this decrease 
of the accretion rate cannot probably change drastically the internal temperature because, at 
high accretion rates, the temperature is close to saturation and depends relatively weakly on the 
accretion rate (see Miralda-Escude et al. 1990 and Fujimoto et al. 1984). Therefore, hopefully, 
our assumption does not change substantially the rate of the magnetic field decay on a spin-up 
line. Note that in our model, the neutron star slides down the corresponding spin-up line with 
the maximal rate but, in reality, the evolution will be slower. The star can leave the spin-up 
line either if accretion is exhausted or if the field becomes too weak to maintain a balance in 
spin up and the field decay. 

3 Numerical results and discussion 

We consider the evolution of the neutron star for a wide range of accretion rates and initial 
magnetic fields. Calculations are performed for the neutron star model with M = 1AM & based 
on the Pandharipande-Smith equation of state (Pandharipande & Smith 1975). This model is 
a representative of stiff equations of state, therefore the radius and thickness of the crust for 
the corresponding neutron star are relatively large, ~16.4 km and ~4.2 km, respectively. We 
choose for calculations one of the stiff models because the magnetic evolution of isolated neutron 
stars based on such models is in a better agreement with the available observational data on 
pulsar magnetic fields (Urpin & Konenkov 1997). The evolution of neutron star models with 
equations of state softer than that of Friedman-Pandharipande seem to be inconsistent with B-P 
distribution of radio pulsars. 

The conductive properties of the crust depend on its chemical composition which, unfortu- 
nately, is rather uncertain for neutron stars. This dependence is not significant during the initial 
evolution when the magnetic field is confined to the layers with not very high density and the 
difference between existing models of the chemical composition is not very large. However, the 
effect of a composition may be important during the phase IV when the neutron star is essentially 
heated and the conductivity is determined by phonon scattering which is strongly sensitive to 
the chemical composition. Therefore, in our calculations, we adopt the simplest model in accor- 
dance with which the crust is composed of nuclei processed in various nuclear transformations in 
accreting neutron stars (see Haensel & Zdunik 1990) during the whole evolution of the neutron 
star. 

In calculations, the impurity parameter, Q, is assumed to be constant throughout the crust 
during the whole evolution and to range from 0.001 to 0.1. The initial spin rotation is assumed 
to be relatively fast, Po = 0.01 s. Note, however, that the results are not sensitive to the value of 
Pq. On the contrary, the evolution is strongly dependent on the original magnetic configuration 
of the neutron star. The initial magnetic field is assumed to be confined to the outer layers of the 
crust with densities p < pq. The calculations presented here are performed for po ranged from 

4 x 10 11 to 10 13 g/cm 3 . These values correspond to the depth from the surface f»820 and 1100 
m, respectively. The choice of po is imposed by comparing the calculated magnetic and spin 
evolution of isolated neutron stars with a crustal magnetic field with the available observational 
data on P and P for radiopulsars (see Urpin Sz Konenkov 1997). The initial field strength at 
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the magnetic pole, B , is taken within the range 10 > B > 10 G. 

The present study is mainly addressed to neutron stars entering binary systems with a 
low-mass companion. For such secondary stars, the main-sequence life-time, t ms , may last as 
long as 10 9 — 10 10 yr. During the main-sequence evolution of the secondary, the accretion rate 
onto the neutron star is determined by both the characteristics of stellar wind and separation 
between the stars, and may generally vary within a wide range. In calculations, we suppose 
M = 10~ 13 , 10~ 15 and lO~ 17 M /yr during the propeller and wind accretion phase; V w = 500 
km/s, a = 5 x 10 12 cm. We also assume £ = 0.1 for the efficiency parameter. Note, however, 
that our calculations show a weak sensitivity of the evolutionary tracks to the particular choice 
of £. After the secondary leaves the main-sequence and fills its Roche-lobe, the accretion rate 
has to be strongly enhanced, M ~ 10 -9 — lO~ lo M0/yr. We do not consider a higher accretion 
rate because of the lack of calculations on the thermal structure of accreting neutron stars with 
M > 3 x 10- 10 M o /yr. 

Figures 1-4 show examples of the magnetic, spin and thermal evolution of the neutron star 
in a binary system for different values of the parameters. Numbers in the top panels indicate 
the evolutionary phase. Note a difference in the horizontal scales for different phases. The 
evolutionary tracks are shown in solid lines if the accretion rates due to Roche-lobe overflow is 
M = lO~ lo M /yr and in dashed lines if M = lO~ 9 M /yr. 

In Fig.l, we plot the evolution for a relatively strong initial magnetic field, Bq = 10 13 G, 
and for po = 10 13 g/cm 3 . The impurity parameter is Q = 0.01, the duration of the main- 
sequence life-time of a low-mass companion is 3 x 10 9 yr. Accretion due to Roche-lobe overflow 
is assumed to last as long as 10 8 yr thus the total amount of accreted mass can reach O.IMq 
if M = lO _9 M0/yr. During the phase III, the accretion rate is lO _15 M0/yr. This model gives 
an example of the neutron star being processed in all evolutionary stages. The phase I, when 
the evolution of the neutron star is not influenced by a companion, lasts a bit less than 10 7 yr. 
During this short stage, the field does not decay significantly thus its surface strength is about 
2 x 10 12 G when the neutron star enters the "propeller" phase. An efficient spin down caused 
by a strong magnetodipole radiation leads to a rapid increase of the spin period from 10 -2 s to 
~ 6 s. During the next evolutionary stage when the neutron star works as a propeller, the field 
continues to decay slowly because, after 10 7 yr of the evolution, the temperature falls down to 
a very low value and the conductivity of crustal matter is high. Note that during this stage the 
conductivity is determined by the impurity scattering and depends on the impurity parameter 
Q. The phase II is relatively long (« 230 Myr) but the field is still very strong (~ 10 12 G) to 
the end of this phase. Due to such a strong field, the neutron star can eject the wind matter 
penetrating into a magnetosphere as long as the spin period is shorter than «4x 10 3 s. When 
a spin rotation slows down to this period, the matter from the wind can fall on to the surface 
of the neutron star and it becomes a soft and not very bright X-ray source. Nuclear burning of 
the accreted material heats the neutron star interior to T ~ 2 x 10 6 K. However, this heating 
is too weak to accelerate appreciably the field decay thus the magnetic evolution during the 
wind accretion phase does not differ essentially from that of an isolated star. Note that the 
wind accretion can accelerate the decay appreciably only if the accretion rate is larger than 
10~ 13 MQ/yr. To the end of the wind accretion, the field strength is still as high as 2 x 10 11 
G. The neutron star can slowly spin up getting some amount of angular momentum from the 
accreted wind matter but this amount is not enough to spin up to a short period, thus the period 
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turns out to be as long as ~ 10 3 s after 3 x 10 9 when the companion ends its main-sequence 
life. The phase III is the most extended for the considered model, however, as it will be shown, 
it may not be the case for other models. The most dramatic changes are experienced by the 
neutron star during the phase IV when an enhanced accretion due to Roche lobe overflow starts. 
Accretion with the accretion rate M = 10~ 10 — lO _9 M0/yr or higher heats the neutron star to 
the temperature T ~ (2 — 3) x 10 s K and decreases substantially the crustal conductivity. If the 
field is reduced only by a factor ~ 50 during the previous 3 x 10 9 years, the accretion-driven 
mechanism leads to a much more efficient decay. Thus, in the case M = lO~ 9 M /yr, the field 
weakens ~ 30 times after 10 8 yr of accreting when the neutron star accretes O.IMq. For a higher 
accretion rate, the decay is even stronger. Since the accreted matter carries out a large amount 
of the angular momentum, the neutron star can accelerate its spin rotation to relatively short 
periods during the phase IV. The spin evolution is non-monotonous (this is shown in the small 
panel in Fig.l). At the beginning of accretion, the neutron star spins up very rapidly and after 
~ 10 4 yr it reaches the corresponding spin-up line. Since the magnetic field is still sufficiently 
strong (~ (1 — 2) x 10 11 G), the corresponding critical period, P eq , is relatively long (~ 1 s). 
The further spin evolution proceeds slowly because the neutron star slides down the spin-up 
line and the rate of this migration is determined by the field decay. Finally, when accretion is 
exhausted after 10 8 yr, the neutron star processed in the all above transformations will work as 
a radiopulsar with a low magnetic field (~ (1 — 3) x 10 10 G) and a short period (~ 0.04 — 0.3 s). 

Fig. 2 shows one more example of the evolution which is typical for stars with a relatively 
weak initial magnetic field. The initial field is assumed to be 10 11 G, other parameters are the 
same as in Fig.l. In this case, the neutron star does not pass throughout all evolutionary stages: 
accretion of matter from the wind is impossible because rotation spins down very slowly due to 
a weak magnetic field and, after 3 x 10 9 yr, the period is still shorter than the corresponding 
value of P eq . A weak initial field results in a very long duration of the phase I for this model: 
the period slows down to a value ~ 0.4 s which allows the infalling matter to interact with the 
magnetosphere only after ~ 600 Myr. However, the magnetic field turns out to be reduced only 
by a factor ~ 10 and is ~ 10 10 G after this long phase. The rest of the main-sequence life-time 
of a companion (which is about 80% of this life-time), the neutron star works as a propeller 
ejecting the wind plasma. Despite of a long duration of this phase, both the magnetic field and 
spin rotation do not change significantly to the beginning of Roche-lobe overflow: the surface 
field strength decreases approximately by a factor 3, whereas the period increases from 0.4 to 
0.46 s. Note a great difference in the periods just before the accretion phase for this model and 
the model plotted in Fig.l where P ~ 10 3 s at t = 3 x 10 9 yr. An enhanced accretion, however, 
can accelerate a spin rotation very rapidly thus an accreting neutron star forgets shortly about 
its period before the accretion phase. Recycling is particularly fast during first ~ (4 — 10) 
Myr of accreting when the neutron star spins up to 10-30 ms depending on the accretion rate 
and approaches the corresponding spin-up line. The further spin evolution on the spin-up line 
proceeds slower but, nevertheless, the star can reach millisecond periods if a mass transfer lasts 
as long as 10 8 yr. For this model, the "recycled" periods when accretion is exhausted are ~ 1 and 
~ 10 ms for M = 10~ 9 and 10~ 10 M Q /yr, respectively. The magnetic field is also substantially 
reduced during the phase IV under the action of the accretion-driven mechanism. To the end 
of this phase, the surface field may be as low as 10 8 G if M = 10 _9 M Q /yr (when the total 
accreted mass is O.1M ) and 3 x 10 8 G if M = 10~ 10 M o /yr (the accreted mass is 0.01M Q ). 
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When accretion exhausts, the neutron star can be observed as a millisecond radio pulsar with 
the period ~ 1 — 10 ms and magnetic field ~ (1 — 3) x 10 8 G. 

In order to illustrate the effect of the impurity content, we plot in Fig. 3 the evolution of the 
neutron star with a strongly polluted crust, Q = 0.1. We also assume that for this model the 
field is initially confined to the layer with a smaller density, po = 10 12 g/cm 3 , corresponding 
to the depth from the surface ~ 900 m. Like the model plotted in Fig.l, this one passes also 
throughout all evolutionary stages. However, the field decays faster for this model because the 
crustal conductivity is lower due to a higher impurity content and the field is initially anchored in 
the layers with a lower density (and, hence, a lower conductivity). A faster decay leads to slower 
spin down and longer duration of two first phases in comparison with Fig.l. Thus, the "isolated 
pulsar" and "propeller" phases last « 13 and 650 Myr, respectively, whereas the duration of 
these phases in Fig.l is ps 6 and 230 Myr. The maximal spin period which is reached to the end 
of the propeller phase does not exceed, however, ~ 300 s. An increase of the impurity parameter, 
Q, causes a particularly rapid field decay during the phases II and III, when the conductivity 
is determined by scattering on impurities. During this two phases, the field is reduced two 
orders of magnitude and is a bit lower than 10 10 G when the companion fills its Roche-lobe 
and an enhanced accretion starts. Accretion of matter during 10 8 yr leads to a formation of a 
millisecond pulsar with the period within the range ~ 2 — 20 ms and with the magnetic field 
(3 — 10) x 10 8 G, depending on the accretion rate during the phase IV. 

Fig.4 represents the evolution of the neutron star with an intermediate initial magnetic field, 
Bq = 10 12 G, and with a longer main-sequence lifetime of the companion, t ms = 5 x 10 9 yr. 
Other parameters are the same as in Fig.l. The evolution does not differ in the main from that 
plotted in Fig.l except an essentially weaker magnetic field of the neutron star at the beginning 
of the accretion phase. For this model, B « 2 x 10 10 G at t = 5 x 10 9 yr. When accretion 
starts, such a weak magnetic field is probably unable to channel the inflowing matter towards 
the magnetic poles, thus the rotation of the neutron star will not cause the observed X-ray 
emission to have a pulsed appearence. Note that the same concerns also the models plotted in 
Fig.2 and Fig.3. On the contrary, the neutron star with the evolutionary scenario represented 
in Fig.l can be observed (at least, at the beginning of the accretion phase) as a pulsating binary 
X-ray sources. The accretion-driven field decay causes a further decrease of the magnetic field 
thus, to the end of the phase IV, it may be as weak as 10 9 G if M = lO _9 M0/yr. If accretion 
lasts 10 8 yr, the neutron star can be easily recycled to a very short period. For this model, the 
final spin periods are ps 10 and 100 ms after accreting 0.1 and O.O1M , respectively. Note that 
the field strength and period of neutron star processed in accordance with the scenarios shown 
in Fig. 2-4 are very similar to those of millisecond pulsars. 

Fig.5-10 show the evolutionary tracks of neutron stars in the B-P plane for different values of 
parameters. We also plot in these figures the observational data on isolated pulsars (dots) and 
pulsars in binary systems (open squares) taken from the catalogue by Taylor, Manchester & Lyne 
(1993). The dotted lines show spin-up lines corresponding to different accretion rates (numbers 
near these lines indicate the logarithm of the accretion rate). The solid lines represent the 
evolutionary tracks. All tracks are calculated for the case of accretion due to Roche-lobe overflow 
with M = lO~ 9 M0/yr. Numbers near the tracks indicate the logarithm of the age, numbers 
with the subscript (a) mark the time required for the neutron star to reach the corresponding 
spin-up line after the end of main-sequence evolution. Large crosses show the points on tracks 
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where a transition takes place between the phases I and II, small crosses with the corresponding 
numbers mark the end of tracks after 10 8 yr of an enhanced accretion. 

Fig. 5 plots the evolution of the neutron stars which accrete matter from the stellar wind with 
different accretion rates. Calculations are performed for M = 1CT 12 (curve 1), 10~ 13 (2), 1CT 14 
(3), 10~ 15 (4) and 10~ 17 Mq/jt (5), the duration of main-sequence lifetime of the companion is 
3 x 10 9 yr. The initial field is Bq = 10 13 G, the density penetrated by the initial field po = 10 12 
g/cm 3 , and the impurity parameter Q = 0.1. The neutron stars starting the evolution with such 
parameters pass typically throughout all evolutionary stages. The only exception is the model 
5 with a very weak stellar wind, M = 10 _17 MQ/yr. This wind is insufficient to slow down 
the spin period to the critical value (11) when accretion from the wind becomes allowed and, 
therefore, the neutron star is still in the "propeller" phase when the companion fills its Roche- 
lobe and an enhanced accretion starts. The duration of phases is obviously strongly dependent 
on the rate of wind accretion. For instance, the initial phase when the neutron star evolves 
like an isolated neutron star lasts as long as ~ 160 Myr if M = lO~ 17 M0/yr and only 10 5 yr 
if M = 10~ 12 M(7)/yr. To the end of this phase, the magnetic field and period are ranged from 
2 x 10 11 to 10 12 G and from 6 to 1 s, respectively, for different models. After this phase, the field 
and spin are reduced sufficiently thus the wind matter can penetrate into the magnetosphere, and 
the neutron star works as a propeller. The "propeller" phase lasts until the neutron star moving 
to the right and down in the B-P plane reaches the corresponding spin-up line. Evidently, the 
star which experiences a stronger accretion from the wind requires a shorter time to approach 
the spin-up line. Thus, one has P = P eq after 5 and 630 Myr of evolution if M = 10~ 12 and 
lO~ 15 M /yr, respectively. As it was mentioned above, in the case of a very weak accretion 
with M = lO" 17 M /yr the spin-up line can never be reached. Accretion from the wind starts 
when the magnetic field is still rather high for the models 1-4, B ~ 10 11 — 10 12 G. This field 
is certainly able to channel the accreting matter towards the magnetic poles thus, probably, 
some of neutron stars accreting matter from the wind in low-mass binaries can be observed as 
weak pulsating X-ray sources. A further evolution of the neutron star on the spin-up lines is 
completely determined by the field decay which can be accelerated due to accretion induced 
heating for a relatively high accretion rate. The wind accretion phase is the most extended 
phase for all models plotted in Fig. 5 with the exception of the model 5. The parameters of 
stars can be drastically changed during this long phase if the accretion rate is sufficiently large. 
For instance, the field is reduced by a factor ~ 10 4 and is only ~ 10 s G for the model 1 with 
the accretion rate M = 10~ 12 M Q /yr. Note, however, that at the end of the main-sequence 
lifetime of the companion the magnetic fields of the models 2-5 which experience accretion with 
M < lO~ 13 M0/yr lie within a relatively narrow range ~ (2 - 10) x 10 9 G despite of a large 
difference in a duration of different phases for these models. The main reason for this is a 
relatively weak heating caused by accretion with such low accretion rates what makes inefficient 
the accretion-driven mechanism at M < lO~ 13 M /yr. 

When the companion fills its Roche-lobe and an enhanced accretion starts, the neutron star 
spins up initially very rapidly: it requires only 0.2 Myr to reach the spin-up line for the model 5 
and 1.6 Myr for the model 2. These time intervals are too short for an appreciable decay of the 
magnetic field, therefore, the corresponding pieces of tracks are represented by approximately 
horizontal plateaus. Note that the magnetic field of the model 1 is very weak after 3 x 10 9 yr of 
evolution (~ 10 8 G), and this model will never reach the spin- up line. Like the case of accretion 
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from the stellar wind, in its further evolution, the neutron star slides down the spin-up line 
and the rate of this sliding is determined by the field decay. However, in the case of Roche- 
lobe overflow, the accretion-driven mechanism is much more efficient and, hence, the decay is 
much faster. During the phase IV, the magnetic field decreases by a factor ~ 30 for the models 
2-5, thus, when accretion is exhausted, all these neutron stars will work as pulsars with weak 
magnetic fields ~ (1 — 5) x 10 8 G and short periods ~ 2 — 10 ms. The model 1 is also spun up 
to a millisecond period but its magnetic field is substantially smaller than 10 8 G. 

In Fig.6, we plot the same as in Fig.5 but for the magnetic field confined initially more deep 
layers, p$ = 10 13 g/cm 3 , and for the crust with a smaller impurity content, Q = 0.01. Evidently, 
the magnetic evolution is much slower for a less polluted crust. Like the case shown in Fig.5, 
the models 1-4 pass throughout all evolutionary phases whereas the model 5 with a very weak 
accretion from the wind (M = 10~ 17 M Q /yr) avoids the phase III. Due to a slower field decay 
and, hence, a faster spin down, the duration of the phases I and II is shorter for the models in 
Fig.6 than for the models with the same wind in Fig.5. Therefore, the wind accretion starts 
when the magnetic field of the neutron star is very strong, B ~ 2 x 10 12 G for all models. As 
in the previous case, the wind accretion is the most extended phase for the models 1-4. To the 
end of the main-sequence evolution of the companion (t ms = 3 x 10 9 yr), the magnetic field of 
neutron stars which experience the wind accretion with a low accretion rate, M < 10 _13 MQ/yr, 
lies in a narrow range from 4 x 10 10 to 3 x 10 11 G. Only accretion with M > 10^ 13 M Q /yr can 
substantially change the magnetic evolution during the phase III. 

Due to a higher magnetic field at t = 3 x 10 9 yr, the models plotted in Fig.6 reach the 
spin-up line on a much shorter time scale when an enhanced accretion starts. Thus, it requires 
only 4 x 10 3 yr for the model 5 and ~ 1 Myr for the model 1. Note that, in the previous case, 
the model 1 could never reach the spin-up line during the accretion phase. If accretion due to 
Roche-lobe overflow lasts 10 8 yr, the models 2-5 which experienced a weak accretion from the 
wind evolve to relatively weakly magnetized pulsars with B ~ (2 — 10) x 10 9 G rotating with 
the period ~ 0.03 — 0.1 s. The model 1 with a very high accretion rate during the phase III 
evolves to the millisecond pulsar with P k, 3 ms and B « 10 8 G. Obviously, the magnetic field 
and period are larger for all pulsars if the accretion phase is shorter. 

Fig. 7 compares the evolution of neutron stars with different initial magnetic fields and with 
the same po and Q. Calculations have been performed for po = 10 13 g/cm 3 , Q = 0.01 and 
tms = 3 x 10 9 yr. It is apparent from this figure that the evolution is strongly sensitive to the 
initial field strength. The weaker is the magnetic field, the longer is the duration of the "isolated 
pulsar" and "propeller" phases and, correspondingly, the shorter is the wind accretion phase. 
For instance, the neutron star can work as an "isolated pulsar" only 0.5 Myr if Bq = 10 G 
but this stage may be as long as 63 Myr if the initial field is weak, Bq = 10 11 G. The strongly 
magnetized star with Bq = 10 13 G can accrete matter from the stellar wind after « 12 Myr 
whereas the weakly magnetized one with Bq = 10 11 G after 600 Myr. To the be ginning of the 
wind accretion, magnetic fields of the considered models lie within a wide range from 10 10 to 
2 x 10 12 G. The corresponding periods are ~ 10 — 1000 s. The wind accretion phase is the most 
extended phase for all models plotted in Fig. 7. During two previous phases the magnetic field 
is typically reduced by a factor ~ 5 — 10, whereas the decrease during a wind accretion is about 
10-50 times. When the companion fills its Roche-lobe and a heavy mass transfer starts, the field 
of the neutron star is « 10 9 and « 4 x 10 10 G for the most weakly and most strongly magnetized 



13 



among the models presented in Fig. 7, respectively. Probably, these fields are too weak to channel 
the accreting matter towards the magnetic poles, thus the neutron star is unable to work as a 
pulsating X-ray source during the phase IV. An enhanced accretion accelerates the field decay 
thus, when the mass transfer is exhausted, the field and period of the models 1-3 turn out to be 
quite suitable for millisecond pulsars. The model 4 which has initially the strong magnetic field 
Bq = 10 13 G may also be transformed into a weakly magnetized and rapidly rotating pulsar to 
the end of the accretion phase. However, the field of this pulsar (~ 2 x 10 9 G) is a bit stronger 
and the period (~ 30 ms) is a bit longer than it is typical for millisecond pulsars. Evidently 
that, if accretion lasts longer than 10 8 yr, the model 4 can also evolve to a millisecond pulsar 
with the standard parameters. 

In Fig. 8, the dependence of evolutionary tracks is shown on the duration of main-sequence 
lifetime of the companion, t ms . Since we address mainly the evolution of neutron stars entering 
binary systems with a low-mass companion, calculations have been done for a sufficiently long 
t ms , 10 10 > t ms > 10 9 yr. Despite of the fact that we plot in this figure tracks for the neutron 
star with a relatively polluted crust with Q = 0.03 and for the initial magnetic field Bq = 3 x 10 12 
G which is certainly weaker than measured fields of many isolated pulsars, it is seen that crustal 
currents can maintain a sufficiently strong magnetic field during the whole evolution. This 
conclusion is in contrast with the widely accepted opinion that neutron star models with the 
crustal magnetic field are not suitable to account for a relatively strong magnetic field of many 
long-lived pulsars. Even for a "polluted" crust, the decay turns out to be very slow thus the 
neutron star can possess the field stronger than 10 9 G after 10 10 yr of evolution and an enhanced 
accretion. 

During the "isolated pulsar" phase which lasts « 20 Myr, the field weakens by a factor « 5 
whereas during the " propeller" phase which is as long as « 630 Myr the field is reduced only by 
a factor 2. Working during « 610 Myr as a propeller, the neutron star spins down to the period 
of the order of 10 3 s which is slow enough to accrete the wind matter. At this moment, the 
field is rather strong, B«3x 10 11 G, thus the accreted matter has to be channeled towards the 
magnetic poles. Perhaps, during the wind accretion phase such a neutron star can be observed 
as a weak and pulsating X-ray source. When the companion fills its Roche-lobe, the field of the 
neutron star can range from 3 x 10 10 to 2 x 10 11 G depending on the duration of main-sequence 
lifetime. The corresponding spin periods lie within the range 100-500 s. Due to a relatively 
strong magnetic field, a spin-up is initially very fast thus the neutron star approaches the spin- 
up line shortly after the beginning of Roche-lobe overflow. The field does not change practically 
during this short time interval (~ 10 4 — 10 5 yr). In the course of the further evolution alongside 
the spin-up line, the field is decreased by a factor ps 30 under the action of the accretion-driven 
mechanism thus to the end of the accretion phase a pulsar can be formed with B ~ (2 — 5) x 10 9 
G and P ~ 20 — 70 ms. If accretion lasts longer than 10 8 yr and the total amount of accreted 
mass is larger than 0.1M@, all the models plotted in Fig. 8 can evolve to "standard" millisecond 
pulsars with B ~ 10 s — 10 9 G and P ~ 1 — 10 ms. Note that if the mass of a low-mass companion 
is relatively large thus the main-sequence lifetime is closer 10 9 yr rather than 10 10 yr, then the 
field at the beginning of the phase IV is sufficiently strong (B ~ 2 x 10 11 G) to channel the 
accreted matter towards the magnetic poles, and the neutron star can work as a pulsating X-ray 
source. To the end of accretion, the field has to be reduced to a very low value insufficient to 
produce X-ray pulses but sufficient for a neutron star to work as a radiopulsar. 
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Fig. 9 shows the dependence of the evolution on the impurity "pollution" of the crust. Calcu- 
lations have been performed for the values of the impurity parameter, Q, within the range from 
0.001 to 0.3. The duration of the main-sequence lifetime of the companion is 3 x 10 9 yr. During 
the initial evolution, the field decay does not practically depend on the impurity parameter. This 
is due to the fact that the neutron star is relatively hot initially and the crustal conductivity 
is determined by phonon scattering which is independent of Q. Later on, when the star cools 
down and the impurity scattering dominates the conductivity, the latter becomes sensitive to a 
"polution" of the crust and the decay may be essentially different for different Q. 

The behaviour of all models plotted in Fig. 9 is more or less standard during the main- 
sequence evolution of the companion. The only exception is the model 1 with extremely large 
Q. Due to a large impurity content and, as a consequence, a low conductivity, the magnetic field 
of this model decays too fast. During the phases I and II which last together « 1.6 x 10 9 yr, the 
field strength has to be reduced more than 1000 times and, when the neutron star approaches 
the spin-up line corresponding to the wind accretion rate M = 10~ 15 MQ/yr, its field is too low 
(B ~ 2 x 10 9 G) to maintain a balance in spin up and the rate of the field decay. Therefore, 
the neutron star does not slide down the spin-up line in its further evolution but moves down 
the B-P plane much faster. To the beginning of Roche lobe overflow the field of the model I 
weakens to a very low value ~ 10 s G. An enhanced accretion causes a further field decay thus 
when accretion is exhausted the field strength of this model is extremely low (~ 10 6 G) and the 
pulsar is practically unobservable. 

The models 2-5 with Q = 0.001 — 0.1 indicate more or less reasonable behaviour. Eventually, 
the neutron star passes throughout all evolutionary phases and, after all transformations, it 
forms a pulsar with a weak magnetic field (B ~ (2 — 30) x 10 8 G) and a short period (P ~ 5 — 40 
ms). Note that, at the end of accretion, the field of the model 5 turns out to be lower than 
that of the models 3 and 4 despite of a smaller impurity parameter for the model 5. Due to 
a substantially higher conductivity, the field of the model 5 diffuses inward very slowly thus, 
after 3 x 10 9 yr, the current maintaining the magnetic configuration is located at much lower 
densities than for the models 3 and 4. When accretion due to Roche-lobe overflow heats the 
neutron star to the temperature ~ 2 x 10 8 K, the conductivity of layers occupied by the current 
becomes essentially lower for the model 5 because of a low density. Therefore, the field decays 
much faster for this model during the accretion phase and, at the end of evolution, it is even 
lower than for the model 3 and 4. 

The dependence of the evolutionary tracks on the initial depth penetrated by the field is 
shown in Fig. 10. We plot the tracks for four values of po : 10 (curve 1), 10 12 (2), 10 13 (3) and 
10 14 g/cm 3 (4). Evidently, the smaller is the initial depth penetrated by the field, the faster 
is the decay. For the model 1 with the magnetic field anchored initially in a relatively narrow 
surface layer with the depth < 600 m, the decay is so rapid that the neutron star cannot even 
be processed in all evolutionary transformations: after 3 x 10 9 yr of evolution the rotation is 
still fast enough to prevent accretion from the stellar wind. Due to this, the model 1 evolves 
from the "propeller" phase directly to the enhanced accretion phase missing the wind accretion. 
However, even for this model, the field at the end of the mass transfer is sufficiently strong thus 
when accretion is exhausted the neutron star can work as a millisecond pulsar with B ~ 3 x 10 8 
G and P ~ 8 ms. The models 2-4 pass throughout all phases and, evidently, the final magnetic 
fields for them are higher and the periods are longer. For example, at the end of accretion, B 
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may be as high as ~ 10 G for the model 4 where the field occupies initially about half a crustal 
thickness. 

4 Conclusion 

We examined the magnetic and spin evolution of neutron stars in close binary systems with a 
low-mass companion. Due to a long lifetime of the companion, evolutionary transformations of 
the neutron star may last as long as 10 9 — 10 10 yr and may be very complex. Likely, most of 
neutron stars in such binary systems experiences a heavy mass transfer ensued due to Roche- 
lobe overflow. This mass transfer can be extremely extended and the total amount of the mass 
accreted by the neutron star can exceed 0.1 — 0.5M Q (see, e.g., van den Heuvel & Bitzaraki 
1995). Besides, the evolution of neutron stars in low-mass binaries has to be influenced by the 
stellar wind during the main-sequence life of the companion. Due to this, the neutron star in 
the course of its evolution passes throughout several essentially different evolutionary phases. 

Calculations presented here are performed in agreement with the standard scenario of evo- 
lution of the neutron star in a close binary system (Pringle & Rees 1972, Illarionov & Sunyaev 
1975). According to this scenario, the neutron star can be processed in four main evolutionary 
phases. Initially, when the star rotates rapidly and the field is strong, the pressure of the pulsar 
radiation does not allow the wind matter to penetrate into the magnetosphere, and the neutron 
star evolves like an isolated pulsar. The duration of this stage may vary depending on the initial 
field strength and the rate of field decay as well as on the intensity of the stellar wind of the 
companion. The initial phase is obviously shorter for a neutron star with a higher magnetic 
field and for a stronger wind. Thus, for the initial field as strong as 10 13 G, the phase I may last 
only ~ 0.1 — 1 Myr if the neutron star captures gravitationally ~ 10 -12 — lO _13 M0/yr from the 
wind plasma. However, this phase may be extremely long (~ 100 — 300 Myr) if the initial field 
is weak, Bq ~ 10 11 G or if the rate of gravitational capture is small, ~ 10 -15 — 10 _17 MQ/yr. 
For the most of considered models, the field is typically reduced by a factor ~ 5 — 10 during this 
phase except the cases when its duration is extremely long. 

When the energy loss due to magnetodipole radiation slows down rotation to such extent 
that the wind matter can interact with the magnetosphere, the neutron star begins to work 
as a propeller ejecting the wind plasma. Obviously, the propeller action is more efficient for 
a stronger magnetic field. The "propeller" phase is usually more extended than the phase I. 
During the " propeller" phase, the magnetic evolution of the neutron star in a binary does not 
differ from that of an isolated star. Since the star is relatively cool to the beginning of the phase 
II, the conductivity of the crustal matter is determined by impurity scattering and may be high 
for a low impurity content. Due to this, typically, the field decays less during a more extended 
"propeller" phase than during a shorter "isolated pulsar" stage. The spin evolution proceeds, 
however, more rapidly because the rate of angular momentum loss is larger for the propeller 
effect than for the magnetodipole braking. During the "propeller" phase the neutron star can 
spin down to a very long period ~ 10 2 — 10 4 s. Note that this period may be shorter for stars 
with low magnetic fields and for a stronger stellar wind. 

In the further evolution, the spin period of the neutron star can reach the critical value 
given by equation (11). When rotation becomes so slow, the centrifugal force is unable to eject 
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the wind plasma penetrating into the magnetosphere, and the matter falls on to the neutron 
star surface. Nuclear burning of the accreted material heats the neutron star, decreases the 
crustal conductivity and accelerates the field decay. Note that an efficiency of the accretion- 
driven mechanism of the field decay is strongly sensitive to the accretion rate. The effect of 
accretion is relatively small for a low accretion rate, M < 10~ 13 M Q /yr, but it may be of a 
great importance for higher accretion rates. During the wind accretion, the neutron star slides 
down the corresponding spin-up line and the rate of this sliding is determined by the rate of the 
field decay. Depending on the duration of this phase and the rate of accretion, the field can be 
reduced by a few orders of magnitude. Since the accreted wind matter carries a some amount 
of the angular momentum, the neutron star can also spin up slightly during this phase but, for 
the most of considered cases, the spin period is still longer 1 s to the end of the main-sequence 
evolution of a companion. Usually, the wind accretion phase is the most extended among all 
evolutionary phases. Note, however, that the star can miss this stage in some cases (low rate 
of mass loss by a secondary star, rapid decay of the magnetic field), thus the evolution goes 
directly from the "propeller" phase to an enhanced accretion. 

In low-mass binaries, the neutron star experiences the most dramatic changes during a 
heavy mass transfer which starts when the companion ends the main-sequence evolution and 
fills its Roche-lobe. Probably, a mass transfer in such systems can last as long as 10 7 — 10 s yr. 
Nuclear burning of the accreted material can heat the neutron star interior to the temperature 
T ~ (2 — 3) x 10 8 K and substantially reduce the crustal conductivity. Due to this, the accretion- 
driven mechanism of the field decay is especially efficient during the phase IV. Typically, the 
field may be reduced about 20-30 times if accretion lasts 10 8 yr, but the decrease may be larger 
for a longer accretion. Depending on the initial magnetic configuration and the duration of the 
main-sequence evolution, the surface magnetic field can fall down to the value ~ 10 8 — 10 9 G 
to the end of a mass transfer. Since the matter accreted from the disc carries a large angular 
momentum, the spin evolves very rapidly. The neutron star reach a spin-up line corresponding 
to the enhanced accretion on a short time scale, ~ 0.01 — 1 Myr. In its further evolution, the 
star moves along the spin-up line while the magnetic field enables to maintain a balance in spin 
up and the rate of the field decay. During this phase, the accretion torque spins up the neutron 
star to a very short period, ~ 1 — 100 ms. 

Our calculations show that the neutron star with the crustal magnetic configuration can 
maintain a relatively strong field during the whole evolution in low-mass binaries. Even if the 
main-sequence life of a companion lasts as long as 10 10 yr and, after that, the neutron star 
experiences a strong accretion with M ~ lO~ 9 M0/yr during 10 8 yr, the magnetic field can be 
sufficiently strong to the end of the accretion phase. When accretion is exhausted, such a neutron 
star works as a radio pulsar with a weak magnetic field, ~ 10 8 — 10 10 G, and a short period, 
P ~ 1 — 100 ms. These parameters are close to those of millisecond pulsars and, perhaps, these 
objects are formed in accordance with the considered scenario. Note that if the main-sequence 
evolution is shorter and a mass transfer is not so extended, the neutron star processed in the 
above evolutionary transformations should have a stronger magnetic field and a longer period 
to the end of accretion. This gives a natural explanation of the origin of those radio pulsars 
in binary systems which lie between millisecond pulsars and the main pulsar population in the 
B-P plane. 
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Figure 1: The magnetic, spin and thermal evolution of the neutron star with Bq = 10 G, 
po = 10 13 g/cm 3 and Q = 0.01. The duration of the main-sequence lifetime of the companion 
is t ms = 3 x 10 9 yr, the rate of accretion from the stellar wind is 10~ 15 MQ/yr. The evolution is 
shown for two values of the accretion rate due to Roche-lobe overflow, M = 10 -10 (solid lines) 
and 10 _9 M Q /yr (dotted lines). The duration of accretion is assumed to be 10 8 yr for both cases. 
The panel shows the spin evolution during the phase IV, t a = t — t ms . 
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Figure 4: Same as Fig.l. but for B = 10 12 G and t rns = 5 x 10 9 yr. 



22 




Figure 5: The evolutionary tracks of the neutron star for different rates of accretion from the 
stellar wind: M = 1(T 12 (curve 1), 1(T 13 (2), 1(T 14 (3), 10~ 15 (4) and l(r 17 M /yr (5). Other 
parameters are: Bo = 10 13 G, po = 10 12 g/cm 3 , Q = 0.1, t ms = 3 x 10 9 yr. The rate of accretion 
due to Roche- lobe overflow is lO~ 9 M /yr, the duration of this phase is 10 8 yr. 
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Figure 7: The evolution of the neutron star with different initial magnetic fields: Bq = 10 
(curve 1), 3 x 10 11 (2), 10 12 (3) and 10 13 G (4). Other parameters are: p = 10 13 g/cm 3 , 
Q = 0.01, t ms = 3 x 10 9 yr. The rate of accretion from the stellar wind and due to Roche-lobe 
overflow is 10~ 13 and lO -9 M /yr, respectively. 
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Figure 8: The evolution of the neutron star in binary systems with a different duration of the 
main-sequence lifetime of the companion: t ms = 10 10 (curve 1), 5 x 10 9 (2), 3 x 10 9 (3), 2 x 10 9 
(4) and 10 9 yr (5). Other parameters are: p = 10 13 g/cm 3 , B = 3 x 10 12 G, Q = 0.03. The 
rate of accretion from the stellar wind and due to Roche-lobe overflow is 10 -15 and lO _9 M /yr, 
respectively. 
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Figure 9: The evolution of the neutron star with different values of the impurity parameter: 
Q = 0.3 (curve 1), 0.1 (2), 0.03 (3), 0.01 (4) and 0.001 (5). Other parameters are: p = 10 13 
g/cm 3 , Bq = 3 x 10 12 G, t ms = 3 x 10 9 yr. The rate of accretion from the stellar wind and due 
to Roche-lobe overflow is 10~ 15 and 10~ 9 M©/yr, respectively. 
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Figure 10: The evolution of the neutron star with the original magnetic field confined to the 
layers of a different density: p = 10 11 (curve 1), 10 12 (2), 10 13 (3) and 10 14 g/cm 3 (4). Other 
parameters are: Bq = 3 x 10 12 G, t ms = 3 x 10 9 yr, Q = 0.01. The accretion rates are the same 
as in Fig. 8 and 9. 
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